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Abstract. - We study the dimensional dependence of the interplay between correlation and 
disorder in two dimension at half filling using 2D t—t' disordered Hubbard model with deterministic 
disorder both at zero and finite temperatures. Inclusion of t' without disorder leads to a metallic 
phase at half filling below a certain critical value of U. Above this critical value U c correlation 
favours antiferromagnetic phase. Since disorder leads to double occupancy over the lower energy 
site, the competition between Hubbard U and disorder leads to the emergence of a metallic phase, 
which can be quantified by the calculation of Kubo conductivity, gap at half-filling , density of 
states, spin order parameter, Inverse participation ratio (IPR) and bandwidth. We have studied 
the effect of disorder on the system in a very novel way through a deterministic disorder which 
follows a Fibonacci sequence. Behaviour of different parameters show interesting features on going 
from a two to quasi one dimensional system. 



Introduction. — Kravchenko et.al. 1 in a series of 
experiments on very weakly disordered 2D semiconduc- 
tors (Si MOSFETS) at very low filling, in the absence of 
magnetic field showed that it is possible for a 2D system 
to undergo Metal Insulator transition. Recently Punnoose 
and Finkelstein have shown that it is possible to identify a 
Quantum Critical point separating a metallic phase stabi- 
lized by electronic interactions from the insulating phase 
where disorder prevails over electronic interactions in cer- 
tain systems [2J. 

Recently the Mott Insulator has been realised for 
fermionic atoms in an optical lattice in 3 dimensions [3] 
whereas most of the previous results were for Bosonic 
atoms. Further the experimental realisation of Anderson 
insulators for Bose gas has been realised [4J by disordering 
the optical lattice site potentials. It is possible to disor- 
der the lattice in a rather controlled manner on optical 
lattices. This has facilitated the controlled experimental 
study of the phase diagram arising due to the competi- 
tion of strong correlations and disorder in these systems, 
causing us to revisit this topic. 

The present work emphasizes on the dimensional depen- 
dence of the interplay between electronic correlation and 
deterministic disorder (and the consequential reflection in 



metal insulator transition(MIT)) as one goes from a two 
dimensional system to a quasi one dimensional system. 
Also the important aspect of the role of temperature in 
this whole scenario has been duly addressed. 

A Hartree Fock mean field treatment of the non- 
disordered t — t' Hubbard Hamiltonian [5] shows Para- 
magnetic to Antiferromagnetic transition at U c = 2.1.e 
The disordered Hubbard hamiltonian was shown to have 
persistent currents [6] , Milovanovic et al [7] investigated it 
numerically and Bhatt et al [5] studied local moment for- 
mation. Dobrosavljevic and co-workers [3] use generalized 
DMFT method to report the interesting fact that increas- 
ing disorder from a clean Mott insulator results in the 
closing of the Mott gap and the stabilization of a metallic 
phase. 

The possibility that a metallic phase may exist in 2D 
due to the interplay of strong correlations and random 
disorder for the half filled disordered Hubbard Model has 
been suggested [10] . 

We solve the problem numerically using Unrestricted 
Hartree Fock(UHF) technique. The UHF method works 
remarkably well even in ID at half-filling and is in ex- 
cellent agreement with results obtained from real space 
renormalization method [H] . In 2D at half filling the UHF 
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method gives result in close proximity with that of quan- 
tum monte carlo method [5]. However, one should mention 
here that the correlated Kondo like processes that lead to 
large effective-mass renormalization in the vicinity of the 
Mott transition can only be captured by a method like 
DMFT and not by the present method [TJ]. The impor- 
tance of inelastic scattering processes for calculating the 
conductivity was shown by a DMFT calculation by Aguiar 
et. aZ[13]. We choose a deterministic binary disorder (that 
follows a Fibonacci sequence), motivated largely by opti- 
cal lattice realisations. Deterministic disorder has been 
one of the forefront areas of research in condensed mat- 
ter physics for the past two and a half decades both in 
experiments [14] and theory [15]. Theoretical studies of 
the quasi-periodic systems have shown that in the non- 
interacting limit the wave function shows a power law lo- 
calization both in one [TB] and two [IT] dimensions. Very 
recent experiment in ID quasi-periodic optical lattice [T5] 
where the system was described by Aubry-Andre Hamil- 
tonian has shown exponentially localized states (Ander- 
son localization) in the large disorder limit. A theoretical 
explanation has been published recently |19) using a Fi- 
bonacci sequence in ID. We benchmark our results and 
trends against old results obtained with random disorder. 
Zero and finite temperature calculations have been car- 
ried out for different values of the disorder strength. The 
impurity concentration remained almost constant as we re- 
duced the dimension, which is an artifact of the Fibonacci 
sequence. We have taken a 40 x 40 lattice which is fairly 
large to rule out finite size effects. 

The system, though insulating in the U > U c regime, 
undergoes a depletion of the charge gap at half filling as 
the disorder is increased, till it reaches a stage where the 
lower energy sites would be doubly occupied in spite of the 
large value of U . A narrow metallic regime emerges while 
going from the Neel state to the state with a distribution 
of doubly-occupied and unoccupied sites. The low tem- 
perature metallic phase gives way to an insulating phase 
upon heating. It has residual spin order(very low value) 
and is thus not a perfect paramagnet. One of our objec- 
tives in the present work is to investigate the evolution of 
the interplay between disorder and interactions and rela- 
tively whose effect gets more dominant as one reduces the 
dimension. The effect of next nearest neighbour hopping 
makes the delocalized phase more robust. 

Disordered Hubbard Hamiltonian. The disor- 
dered 2D t-t' Hubbard Hamiltonian considered by us is as 
follows: 

H = WH a - t J^(c+ c l+Sia + h.c) 

i i 

+*'X)(4 c *+^ + M (i) 

i 

Hint = U ^ n i\ n ii ( 2 ) 
i 

H = H D + H inc (3) 



Table 1: Number of A & B sites for different system sizes. 



size 


A 


B 


A/B 


40 x 40 


989 


611 


1.608 


40 x 14 


346 


214 


1.616 


40 x 10 


247 


153 


1.614 


40 x 6 


148 


92 


1.608 



Here t and t' are nearest and next nearest neighbour hop- 
ping terms, U is the on site Hubbard interaction, ti is the 
site potential at i-th site, while cf a , Ci a , , are cre- 
ation, annihilation operators for the electron of spin cr, and 
the number operators for up and down spins respectively 
at site i. 

Deterministic Disorder. — Deterministic disorder 
is neither periodic nor fully random and have long range 
correlations(correlated disorder). The site potential e; is 
deterministically disordered. It follows a Fibonacci se- 
quence which is generated as A — > AB and B — !> A, 
where A and B are the two different sites. In our case 
site B has a higher site potential than site A (eg > (a, 
W = e B — ea, being the disorder strength). A typi- 
cal Fibonacci chain in a particular generation looks like: 

ABA ABABA ABA ABABA ABABA We have used 

the Fibonacci sequence and generalised the idea of quasi- 
periodicity in one dimensions to two dimension by essen- 
tially wrapping up the one dimensional Fibonacci sequence 
over the two dimensional square lattice. Moreover the 
number of lattice sites need not be a Fibonacci number. 
In our case the N lattice sites picks up the first N entries 
of the sequence from the left, thus allocating the site po- 
tentials. This sort of disorder can also be referred to as 
correlated disorder. 

Discussion of Numerical Method. We solve 
the Hartree Fock decoupled Hamiltonian with an initial 
guess(seed) of and n^, and solve it self-consistently 
till the solutions of and converge to a difference 
of less than 10~ 7 for every site. We settled up with that 
initial seed which minimises the ground state energy. This 
strict convergence criterion has been employed to ensure 
that we stay away from local minimas while exploring the 
energy landscape. This is extremely important as we ap- 
proach the disorder induced localized phase using a vari- 
ational approach like Hartree Fock to treat interactions 
as our only guiding principle in choosing the ground state 
configuration is energetics. In the low and medium dis- 
order regime if we start from the Neel order state as the 
starting seed, we get the desired convergence very quickly, 
while in the high disorder regime near the MIT, we have 
found the desired results starting with a near paramag- 
netic configuration as the initial seed(PM seed). 

Unrestricted Hartree Fock . In UHF theory the Hub- 
bard interaction term is decoupled keeping in the site de- 
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pendence and thus one obtains modified site potentials for 
up and down spins which reads as: 

t i]r = Ci + U < n it >; e it = e< + U < > (4) 

This gives us the following Hamiltonian which is decoupled 
into up and down spin parts: 

i i 

+t'J2( C tc l+ 5 2 a+h.c) (5) 

i 

Where e ia are now defined by eq. (J4j) . 

Quantities calculated. — Once the energy spectrum 
is obtained the energy gap at half-filling is calculated by 
the difference between the energy level at half-filling and 
the one just above it. The expression reads as: 

Egap = E N / 2 +i — En/2 (6) 

where N is the number of sites in the 2D lattice. The 
ground state energy E g is calculated by summing up the 
single particle energy states for both spin species up to the 
Fermi level corresponding to the desired filling. The Fermi 
level at half filling is calculated as the mean of the highest 
filled state and the lowest unoccupied one to ensure that 
the system remains half filled. All finite temperature cal- 
culations are performed with respect to this Fermi energy 
with the appropriate Fermi distribution factors. If and 
Ujj, are converged values of of the occupation numbers of 
up and down spins respectively at the i-th site, then the 
spin order is defined as: 

5 = ^(E (-i)>*t-«4)+ E (-i) (i+1) («,t-^)) 

i,evenrow i,oddrow 

(7) 

For disordered non interacting systems, the Kubo formula, 
at any temperature is given by: 

a H = 4E("«-^)^^-(^- e «)) ( g ) 

with A — ire 2 /ha,Q 1 ao being the lattice spacing, and n a = 
Fermi function with energy e a — /i. The / Qj a are matrix el- 
ements of the current operator j x = it a (c i+xa a Ci, c — 
h.c), between exact single particle eigenstates \tp a ), \ipp), 
etc, and e a , ep are the corresponding eigenvalues. In this 
paper, conductivity/conductance is expressed in units of 
A = Tre 2 /hao- 

We calculate the 'average' conductivity over 3-4 small 
frequency intervals Aw(Au; = nu) r , n = 1,2,3,4), and then 
differentiate the integrated conductivity to get a(u)) at ui = 
nuj r , n = 1,2,3. These three values are then extrapolated 
to uj = to get Ode- For the sake of simplicity, and keeping 
all other consistency checks in mind [20], we have taken 
the value of a(oj r ) as the actual value of aa c - 



The average of c(w) over the interval [0, cu r ] is 

1 f" r 

a in t(uJr, n,N) = — / a(ui, /j,,N)duj (9) 

w r Jo 

io r is set to be sufficiently larger than the average spac- 
ing between energy levels, given by D/N, where D is the 
full bandwidth for the interacting 2D disordered Hubbard 
Model that we are considering and is simply calculated as 
D = En — Ei, where En and E\ refer to the topmost and 
lowermost eigenvalues. For N = 1600, u) r — 0.05. 

In the discussion below, all length scales are normal- 
ized by the lattice parameter a, and we use dimensionless 
energy parameters U, W and t', scaled by the hopping 
amplitude t. 

IPR which gives the measure of how many sites are 
participating in a particular eigenfunction is defined as 
following [21] : 

IPR{E n ) = j2\t ] \VCt\^ l) \ 2 ) 2 (10) 
(=i i=i 

where I corresponds to the site index, E n is the nth eigen- 
value. The cited ref calculated IPR for a non interact- 
ing system with Fibonacci modulated site potential. It is 
clearly visible from the above expression that larger value 
of IPR corresponds to localized electronic states whereas 
small value of IPR indicates an extended state. In our 
case the topmost eigenfunction near the Fermi level was 
considered for this purpose. 

Numerical Results. — We have calculated the gap 
at half filling, charge and spin order, bandwidth, IPR and 
the <t|* of the 2D t — t' disordered Hubbard model at half 
filling as we vary the t', Disorder potential W and also 
reduce the dimensionality along the y direction. We study 
the problem for intermediate disorder strength (W < U) 
where the underlying Neel magnetic ordering(due to U) 
is not killed and is predominantly present, for a 40 x 40 
lattice. We have shown here our calculations for system 
sizes 40 x 40, 40 x 14, 40 x 10 and 40 x 6. We have tabulated 
the number of A and B sites and the ratio A/B for the 
different sizes which is given in table [1] 

In Fig. QJi we study effect of size reduction along y axis. 
We show that for t' = 0.2, E gap is zero till a critical U c = 
2.1, above which it becomes non-zero and matches well 
with Hirsch's result [5] . U c is unchanged all the way down 
to 40 x 14 below which U c starts decreasing slowly at first, 
and below 40 x 8, it shows a steep fall, as the quasi ID 
effect sets in and electronic motion gets constrained. In 
Fig. [Tp Eg decreases almost linearly as we increase t', in 
the regime ( U > U C ,U C = 2.1), as t' kills the Neel ordering 
due to U. 

Fig. [2k shows the plot of the gap (as we reduce n y ) 
for t' = 0.2 against disorder strength W for U = 3.0, 
which renders the undisordered phase antiferromagnetic. 
The disorder reduces the gap as now the electrons of both 
spin species will try to avoid the site with higher energy, 
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Fig. 1: Fig. [IJl shows the effect of reducing n y on the gap at half 
filling. For t' = 0.2, gap is zero till U c = 2.1. U c is unchanged 
all the way down to 40 x 14 below which U c starts decreasing 
slowly at first, and below 40 x 8, it shows a steep fall, as the 
quasi ID effect sets in and electronic motion gets constrained. 
Fig. Hp shows the fall in gap at half filling with increasing il . 
The plot is for W = 0, U = 3 and N = 40 x 40 



thereby allowing some double occupancies on the lower 
energy site. The gap decreases rather slowly for low values 
of W and finally falls sharply at a certain W — W c , for a 
fixed system size. 

We find that for a fixed disorder, though E g increases as 
we reduce n y , the rate of fall in E g increases with increase 
in disorder as we reduce n y , highlighting that effect of 
disorder is more pronounced than the effect of Hubbard 
U as one reduces n y . We will see the effect of this in our 
Kubo conductivity results. 

The corresponding spin order plot also follows a simi- 
lar trend (Fig. Wp)- However the reduction in this case is 
more gradual than the gap at half filling. This is because 
the system prefers to retain the antiferromagnetic config- 
uration to gain maximum kinetic energy due to hopping, 
till a critical value of W upto which U is effective, and 
beyond which the spin order reduces very rapidly. We 
find this critical W ~ 1.75, at which point the Hubbard 
gap roughly diminishes from first to second place of dec- 
imal. The reduction in magnetic ordering is much slower 
than the gap at half filling, and there is significant residual 
magnetic ordering^ ~ 0.25) when A has already gone to 
zero. 

Fig. [5^ shows the plot of the dc conductivity (erg) against 
disorder strength W at T = for U = 3,3.2 and 3.4, which 
is well in the antiferromagnetic regime even for t' = 0.2. 
Due to interplay between correlation and disorder, a a in- 
creases from zero at a critical value of W = W c , whose 
value increases with U. Thus it takes higher values of dis- 
order to form doubly occupied sites as one increases U at 
half filling. This plot cannot be extrapolated indefinitely 
on the lower side as U has to be greater than U c . 

Fig. [3b shows the plot of Kubo conductivity at zero 
temperature as we change t' for the disordered problem. 
It is found that as we increase t' from zero it takes a lower 
value of disorder to first generate a non zero value of the 
Kubo conductivity and then render the system metallic. 

Fig. [3J: shows the plot of the T = <7o against W as 
we reduce the system size along y direction, thus going 
over to a quasi-lD system. We need lower disorder value 




Fig. 2: Fig. [2^l shows the plots of the gap at half filling and 
Fig. [2p shows the plots of spin order for t' = 0.2, U = 3.0 as a 
function of disorder strength W, for various values of n y . 



to annul the effect of U, as we lower the dimensionality 
again highlighting the fact that the effect of disorder, for 
a fixed value of U becomes more robust as we decrease 
the dimensionality from 2D to quasi ID. The faster rate 
of depletion of the Hubbard gap A as we reduce the di- 
mensionality as seen in Fig. is thus connected to the 
o"o — W curves for different n y in this way. 

Fig. show the plots of er against temperature for 
U = 3.2 and shows an insulator to metal transition with 
increasing W at low temperatures. Our calculation shows 
a finite conductivity at zero temperature, which is con- 
sistent with experimental results. The region kT < 0.02 
corresponds to the diffusive sector while for kT > 0.05 we 
are in the ballistic regime. W c corresponds to the value 
of W for which ao at T = becomes non-zero and slowly 
starts to increase, while at W = W m the system becomes 
a metal as the Hubbard gap is completely killed by the 
disorder. In between there is a narrow range of W within 
which ao starts increasing from zero, but the system is still 
however not a metal. Further increasing W, we enter a 
metallic regime where the rate of fall in conductivity with 
temperature decreases. This indicates that it is a highly 
disordered metallic phase (dirty metal) . The ground state 
of the highly disordered insulator is extremely difficult to 
obtain within Hartree Fock, which takes a very long time 
to converge to the required precision. Even when it does, 
there is a possibility of the system getting stuck in one of 
the many local minimas and not being able to access the 
true ground state. We therefore venture upto the point 
where there are no convergence problems (system is metal- 
lic). 

Though the effect of U becomes stronger, with the Hub- 
bard gap becoming larger as we decrease the dimensional- 
ity of the system, it is clear that disorder effect surpasses 
the effect of U as one lowers the dimension as evident 
from Fig. (3Jd,e,f). Fig. [3Jl,e,f shows the finite tempera- 
ture plots of <To for different parametric values of W for 
n y = 40, 14 and 6 respectively. We again observe the 
general behaviour that the dc conductivity of an insula- 
tor/metal rises/falls with increasing temperature in the 
low temperature regime. However the onset of the metal- 
lic phase occurs for lower values of W as we reduce the 
dimensionality. 

We see that as we heat the system, the system goes from 
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Fig. 3: Fig. [3K,b,c are plots of Kubo conductivity against dis- 
order. Fig. [3JI,e,f are plots of Kubo conductivity against tem- 
perature. In Fig. [3^l plots are for U = 3,3.2,3.4 from left to 
right. The Kubo conductivity becomes non-zero and system 
becomes metallic for higher values of W as U increases. The 
curves in Fig. |3p are for t' = 0.2, 0.15, 0.1, 0.05, from left 
to right. Fig. [3ji shows the effect of reducing n y on the Kubo 
conductivity from left to right, n y = 6,10,14,30,40. Fig. [3JI 
shows MIT at low temp driven by increasing disorder. The 
curves correspond to W — 1.8,2,2.1,2.15,2.2 from bottom to 
top. Fig. [3£ and Fig. [3J are for n y — 14 and 6 respectively and 
highlight how the MIT takes place at lower value of disorder as 
we reduce n y . The curves in Fig. [3^ are for W = 1.45 to 1.95 
in steps of 0.1 from bottom to top, while Fig. [3J are for W = 
1.4, 1.85 in steps of 0.05 from bottom to top. 



a low temperature metallic phase to a high temperature 
insulator phase. This effect is further seen to go away as 
we reduce the dimensionality. 

Fig. @Ji,b shows the plot of the density of states(DOS) 
for a fixed value of U = 3 for n y = 40 and 6 respectively, 
highlighting how disorder closes the gap when we reach 
W m . However the gap closes for lower value of W as we 
decrease n y again emphasising the increased effectiveness 
of disorder as we reduce the dimensionality. We compare 
our result with the work done in the first reference of [5]. 
We observe from their LDOS plot that the MIT for the 
random disordered 2D Hubbard model occurs roughly at 
W = 2.7 for U = 1.25. In our work for the 2D t - 1' (for 
t'=0.2), deterministically disordered Hubbard model the 
MIT occurs at W = 1.85 for U = 3 (Fig. gji). Fig. gt,d 
shows the IPR and the bandwidth of the system for a fixed 
value of U > U c , as we increase the disorder strength. 
Fig. 2Ji clearly shows that for the interacting case the ef- 
fect of disorder on band- width is quite non-monotonic and 
non-trivial. The bandwidth increases for very weak(low) 
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Fig. 4: Plots of DOS in Fig. UK and b, IPR in Fig. |IJj and 
BandWidth in Fig. [4ji. The gap closes for W = 1.85 for n y — 
40(Fig. UK) whereas for n y = 6(Fig. |Ud) it closes for W = 1.6. 
The IPR and Bandwidth figs are for t' = 0, 0.1, 0.15 and 0.2 
respectively from top to bottom. 



disorder, then it stops growing and surprisingly starts de- 
creasing with increasing disorder. This is the regime where 
the states at the Fermi energy have already become delo- 
calized as can be seen from Fig. 0t(the IPR curve also be- 
haves non-monotonically). In this regime more and more 
states are introduced in the gap region. Since total num- 
ber of levels are fixed, the bandwidth decreases. 

The localization in the Hartree Fock single particle ex- 
citations for t' = 0, as reflected in the IPR peak below 
W = W m is rapidly eroded as we increase t'. This is seen 
in the rapidly falling peak structure in the IPR curve as we 
increase t' . Thus with increasing t', there is a subtle tran- 
sition in the nature of the wave-functions itself(localized 
to delocalized) as the metallic phase sets in. 

Fig 5(a) represents the plot of gap at half filling against 
disorder, for three different realizations of the Fibonacci 
sequence. One can see from these plots that the results 
remain almost invariant as we change the realization. The 
results shown here are for U = 3 and t' — 0.2 for a 40 x 40 
system. 

Fig 5(b) shows our results of the dc conductivity at zero 
temperature as a function of disorder, in the range where 
the dc conductivity just starts to pick up from zero. We 
have shown the results for 2 system sizes, namely 32 x 32 
and 40 x 40. It is clearly observed that the critical value 
of disorder W c where the dc conductivity picks up from 
zero is almost the same. This shows that the finite size 
effects are already quite small by the time we reach 32 x 32 
system size. 

We believe that our simulation results provides useful 
and direct insight into the competition between correla- 
tions and disorder, which has been introduced in a sim- 
plified(deterministic) manner and is capable of capturing 
all the subtleties involved. It gives a good estimate of the 
Kubo conductivity as one goes from 2D to quasi ID sys- 
tems in presence of disorder and correlation and indicates 
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Fig. 5: Fig. [5^, shows the effect of disorder on the gap at half 
filling, for three different realizations of the Fibonacci sequence, 
for a 40 x 40 system. The uppermost, middle and lowermost 
curves correspond to the Fibonacci sequence beginning with 
first, seventh and third entries of the sequence respectively. 
Fig. [5]d shows the plot of dc conductivity at T — for 32 x 
32(lower curve) and 40 x 40(upper curve) system size. 



the emergence of a metallic phase. As we tune the disor- 
der to a very high value, the system becomes a metallic 
glass (highly disordered metal). 

* * * 
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